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Abstract 


A fuel cell powered system is regarded as a high current and low voltage source. To boost the output voltage of a fuel cell, a DC/DC converter is 
employed. Since these two systems show different dynamics, they need to be coordinated to meet the demand of a load. This paper proposes models 
for the two systems with associated controls, which take into account a PEM fuel cell stack with air supply and thermal systems, and a PWM 
DC/DC converter. The integrated simulation facilitates optimization of the power control strategy, and analyses of interrelated effects between the 
electric load and the temperature of cell components. In addition, the results show that the proposed power control can coordinate the two sources 


with improved dynamics and efficiency at a given dynamic load. 
© 2006 Elsevier B.V. All rights reserved. 
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1. Introduction 


A PEM fuel cell powered system is composed of a PEM 
fuel cell system, a DC/DC converter with a capacitor, and a 
voltage-impressed inverter that controls a propulsion motor. A 
single cell is constructed by stacking different layers composed 
of bipolar plates, gas diffusion layers, catalyst and membrane. 
When hydrogen and oxygen chemically reacts in the catalyst, 
DC power is generated with heat and water as byproducts. Fig. 1 
shows a typical fuel cell powered system with a PWM DC/DC 
converter, storage and an inverter driving the propulsion motor 
for mobile applications. 

The open circuit voltage of a single cell is theoretically 1.23 V, 
but drops when a load is applied. The rapid output voltage drop is 
caused by increased over-potentials at the electrodes and mem- 
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brane. Thus, the fuel cell exhibits a typical characteristic of alow 
voltage and a high current power source. If a fuel cell system 
should drive an electric motor via an inverter, the output voltage 
of the fuel cell should be stepped up to a level of the DC bus volt- 
age that is at least equal to the peak value of the output voltage of 
the inverter. Boosting the voltage can be electronically carried 
out by a DC/DC PWM converter that, in principle, stores energy 
in an inductor when a switch is turned on, and discharges the 
energy in the inductor to a capacitor when the switch is turned 
off. Once the switching period is properly controlled, the voltage 
of the fuel cell can be boosted to the voltage level at the capaci- 
tor that can be kept constant at varying loads [1-3]. Therefore, a 
DC/DC converter is an indispensable part of the power system 
to ensure a power flow from the fuel cell to the electric machine 
for securing the performance and high efficiency. 

The balance-of-plant (BOP) is the ancillary part that supplies 
fuels and coolant, also and removes the byproducts. The BOP can 
be divided into four subsystems: for hydrogen, air supply, water 
supply and heat removal. The fuel delivery system consists of a 
tank for hydrogen storage, while the oxygen is delivered by air 
through a compressor or a blower with humidifier to maintain the 
humidity in the cell. The heat generated is removed by a coolant 
circuit that includes a radiator, pump, condenser and a reservoir. 
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Nomenclature 

A area (m2) 

C mass concentration (kg m7?) 
C capacitor (F) 


Cp specific heat (Jkg~! K7!) 

F Faraday number 

Fr fontal area (m?) 

h heat transfer coefficient (W m~? °C7!) 
i current density or current (A cm~? or A) 
m mass (kg) 

M molar mass (kg mol`!) 
n number of cells 

N mole flux (mol s~! m7?) 
Nol number of poles 


power (W) 
Q heat transfer (J) 
Relec electrical resistance (Q) 
s entropy (Jmol~! K7!) 
s balancing factor 
t thickness (m) 
v, V voltages (V) 
W mass flux (kg s7! m7?) 


Superscripts, subscripts 


an anode 

bl blower 

ca cathode 

cv control volume 
8 gas 

i index 

membr membrane layer 
S stator 


Greek letters 


E porosity 

T tortuosity 

n efficiency 

0) flux linkage (V s rad™!) 


In addition, a load management is necessary to coordinate all 
aforementioned components to deliver the load power required. 
Four different control systems are therefore needed to obtain 
reliable operations of the fuel cell system as a power source. 

This paper will describe the models for the fuel cell stack, 
air supply system, thermal circuit and a DC/DC PWM con- 
verter as key components in the system. The associated controls 
are briefly introduced. Finally, a new power control strategy is 
proposed. 


2. Modeling of a fuel cell stack, air supply and thermal 
circuit 


Fig. 2 shows a block diagram of a fuel cell powered system 
that includes a stack and four subsystem: hydrogen supply, air 
supply, water and thermal managements. 


If it is assumed that the reactants are fully humidified and 
the hydrogen flow rate supplied is proportional to the air flow 
rate, the power of the stack can be only controlled by a flow rate 
of the air that is supplied by an electrically driven air blower. 
Heat produced by the stack is extracted by the coolant that in 
turn exchanges it with a radiator. Thus, the model for a fuel cell 
powered can be reduced to a stack with an air blower and the 
coolant circuit, and the DC/DC converter. 


2.1. PEM fuel cell stack 


Most fuel cell models, which describe physical behaviors of 
a PEM fuel cell, are described either by empirical equations 
fitted to the curve of the polarization characteristics or by com- 
putational fluid dynamics (CFD) to solve the mass and charges 
transport. The former can be used to describe a steady state 
behavior, but lacks the full dynamics, which is complemented 
by reflecting the double layer capacitor [4]. However, these mod- 
els do not provide the gas dynamics dominant along flow paths. 
For instance, the porosity in the gas diffusion layer (GDL) causes 
a partial pressure drop that increases the catalyst over-potential 
according to the Butler- Volmer equation. 

On the other hand, the CFD based models have been widely 
employed to analyze physical phenomena in a single cell, but 
limited to represent the dynamic characteristic in conjunction 
with components of the BOP that can be used for analysis of a 
power system. 

Consequently, a new transient model is required. The model 
proposed in this paper is based on empirical equations and 
takes into account three major effects: the water balance in the 
membrane, the gas dynamics in the GDLs and the temperature 
distribution in a cell that is briefly described in the following 
sections. 


2.1.1. Model improvement for the stack 

A cell is constructed by individual model for layers, where 
the ohmic over-potential in the membrane, the activation 
over-potential in the catalyst on the cathode side, and the concen- 
tration over-potential are considered to obtain a current-voltage 
static relationship at a cell. The voltages drops for a cell are 
dependent on the reactant partial pressure, temperature and 
membrane water content. 


veel = E (p, T) — Vact (p, T, i) — Vohm(i, Amembr, T) 
—Veconc(P, T, i) (1) 


Ustack = A X Veell 


The dynamics of the fuel cell system are described with the 
changes of mass flows of the air and water in the fuel cell and 
air supply system. The results show a typical J-V characteristic 
of a single cell with a transient behavior mainly representing the 
mass transport in the air delivery system [5-8]. The dynamic 
behavior can be further improved by three effects: (1) water 
dynamics in the membrane, (2) partial pressure drop in the gas 
diffusion layer (GDL) and (3) temperature variation. 

Water content in the membrane determines the proton con- 
ductivity. The dynamics of the water content is described by 
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Fig. 1. Typical configuration of a PEM fuel cell powered system for automotive applications. 


two effects, the electro-osmotic driving force by the different 
electrochemical potential at the anode and cathode, and the dif- 
fusion caused by the water concentration gradient at the two 
boundaries. Considering the water mass flows at the boundaries 
of the membrane layer, the relationship is expanded as follows 
[9-11]: 
CH,O,mass/ MH,0 
(Pdry,membr/Mmembr) = bCy,0,mass/My50 
d(CH,0,mass Acellfmembr) 
dt 
—Wele,membr,ca F Waiff, membr,an =f Waitt, membr,ca 


Amembr = 


(2) 


= Wele, membr,an 


Mvwater,membr 


On the other hand, the reactant entering the cell diffuses 
through the GDL before reaching the catalyst layer that sig- 
nificantly affects the overall dynamics of the reactants. This 
diffusion effect is reflected by using the mass continuity (3) 


Humidifier 
Inlet Manifold 


Blower 
Inlet Manifold 


and the Stefan—Maxwell Eq. (4) [4]: 


Eg OPi ON; 
ZE =0 3 
RT ot dy @) 
3 
= (piNk — pk Ni) (4) 
T? dy DaD ý ai 


Hence, i, ke (1, 3), where pı is the oxygen partial pressure, 
and p2 =Psat(T) and p3 are the water vapor and the nitrogen par- 
tial pressure, respectively. The diffusion coefficients of peaDik, 
Dik,orin and Dix orin(T) include the cathode pressure of pea, sum- 
ming the species partial pressures. The parameter Tt is a constant 
describing the pore curvature of GDL. 

If a cell is assembled with cubical layers, where the thermo- 
physical properties are isotropic and constant, then according 
to the energy conservation equation, the total energy changes in 


SSeS SS ey 


Fig. 2. A schematic diagram of a PEM fuel cell system. 
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Table 1 
Simulation parameters 
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Fuel cell 
n 
Arc (cm?) 


Proton conducting model 
bii 
Dio 
b2 
na 
Dw 


Gas transport model 
Dest (m? s7!) 
Psat 


Electrochemical reaction model 
po (bar) 
Tref (K) 
Erref (V) 
Acatl, eff/Acell 


Thermal model 
H, gas 
CPgas 
Pgas 
Frarea (m°) 
Mres (kg) 
hA (IK7!) 


Thickness (m) Density (W mK~!) 


381 
280 


0.5139 
0.326 

350 
KCwater) 
AT, Cwater) 


AP, T) 
AD) 


1.0 
353.15 
1.229 

SU, T, Poa) 


KP, T) 
SP, T) 
SP, T) 
1 
5 
16.66 


Heat conductivity (Ikg~! K7!) 


[7] 
[7] 
[7] 
[7] 
[7] 


[5] 
[5] 


[19] 
[19] 


[13] 
[13] 
[13] 


Specific heat (kg m~’) 


Geometrical data for layers [13] 


Coolant channel 0.002 1400 
Plate 0.001 1400 
Gas channel 0.001 1400 
GDL 0.0004 2000 
Catalyst layer 0.000065 387 
Membrane layer 0.0001275 1967 


30 935 
52 935 
52 935 
65 840 
0.2 770 
0.21 1100 


a controlled volume are equal to the sum of energy exchange 
at boundaries and internal energy resources. In fact, the energy 
exchanges at boundaries occur by three factors: (a) the mass flow 
into each volume; (b) the conduction heat transfer across the 
cell; (c) the convection heat transfer occurring between bipolar 
plates with the coolant and the reactants. Thus, the thermal- 
dynamic behavior can be described with the following energy 
conservation Eq. (5) [11]: 


SOCP; Pi mass Ace dTev 
- P; Pi, mass Acellfev FP 


E 5 Win AcentCp j(Tin — Tey) + 


Qconv Acell 
< 


fase Howl convection heat transfer 


+ Ò cond A cell 
< 


convection heat transfer 


+ Osou (5) 
WH 


sources 


On the other hand, the internal energy source is composed 
of entropy loss and the chemical energy required for protons to 
overcome the barrier of the over-potentials in both catalyst layers 
(6). In addition, others are ohmic losses caused by a transport of 


electrons and protons in the cell [12]: 


; ; TAs i 
Osou = i Acel TF + n + iAcen Relec (6) 


All models were coded by blocks given in MATLAB/Simulink. 


2.1.2. Parameters and simulation 

The parameters and reference data for the models chosen are 
as follows (Table 1), which are partially empirical [4,6,12,18]. 
In addition, the maximum voltage and current used for this study 
are 0.95 V and 0.95 Acm~?, respectively. 

Multi-run simulations have been conducted to investigate the 
static and dynamic behavior of a single cell. The static behavior 
is analyzed by calculating the typical polarization at different 
temperatures shown in Fig. 3 and temperature distribution at 
a constant current. Only the J-V characteristic is taken at the 
varied temperatures from 333 to 353 K with a step of 10K. As 
the temperature rises, the water removal is eased. Specially, the 
effects are considerably high at the range of the high cell current 
where more water is produced. This result is comparable with 
the CFD analysis [13]. 

Fig. 4 shows the temperature distribution through the cell 
at seven minutes. The temperature changes dynamically and 
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Fig. 3. I-V curve for different cell working temperature cell: p = 1.0 bar. 


consequently the voltage of a cell is no loner assumed to be 
constant at a given load. Particularly, the temperature in the cat- 
alyst on the cathode side shows the highest peak among others 
because of the losses associated with the higher over-potential 


i wp (—20.581(p*)* — 1.4415 x 1073 p* + 4.1333 x 107), 
bl = 


otherwise, 


than the anodic one. The maximum temperature difference 
between the catalyst at the cathode side and the coolant channel 
is 7K. 


2.2. Air supply system 


The air supply system consists of a blower, a humidifier, 
plumbing for the inlet and the outlet. The outlet of a blower 
driven by an electric motor is connected to a humidifier and an 
inlet of flow channels through pipes. In this study, the humidifier 
is simplified as an ideal one without any associated dynamics 
and energy losses. 

The blower is constructed with an impeller driven by an 
electric motor. The dynamic response of the blower is described 
by the inertia of the motor and the impeller, and the torque 
produced by the motor. Hence, the torque produced by the 
motor, Tpi.m, is a function of the stator resistance, Rs bm (Ohm), 
magnetic flux constant, ®pim (V srad~!), and the number of 
the poles, Npim,p1 with the stator voltage, Vim (V), (8) [14]. 
Therefore, the flow rate (W,)) can be controlled by the voltage 


aiaee] 
l-9 T 
3 S| E & || 2) 
Anode 2 gi gia Ei Cathode 
x 
FEAHHK 
ž EJ ajs 
T 354) - 2 
x 
z 353 
& 352) {ltt} td 
E 351 ee ee a >> 
8. 350 e eet Nee 
E 349) {ill tly 
348 _ 


Fig. 4. Temperature distribution in the cell at 7 min after a startup (left to right: 
anode coolant channel to cathode coolant channel). 


of the motor, Vpim (V). 


dap1 1 ( Woi APoinol ) 
= Tbl,m 
dt Joi bl Pambbl 


ee 3 (=n) (==) 
= ee 2 Rs,bl,m 
Noi,m,pl 
x ce = (=en) Pomon (8) 


where wp) is the angular velocity (rad sy, Jp, the rotational 
inertia (kg m2), Nblm the motor efficiency, pp) the pressure (Pa) 
and Pamb is the density (kg im 

The parameters for the blower are derived by characteris- 
tic data and specifications given by Phoenix Analysis & Design 
Technologies (PADT) [15], which includes both the flow param- 
eter and overall efficiency versus the head parameter. The 
equations describing the behavior of the impeller are given as 
follows [14]: 


(7) 


p* <9 x 1074 Pa s? rad~? 


* -3 (9) 
wpi(—1.7973 p* + 1.6409 x 107°) 
Wa \> Wo \7 
Nb) = —2.8831 x 10837) +9.5115 x 108 (=) 
bl bl 
4 { Wo 
+ 1.3087 x 104 ( —> ) + 0.17945 (10) 
bl 


e. 2 
where p` is ((Pea — Pamb)/@p) . 
The inlet and outlet manifold pressures are governed by the 
mass conservation equation: 


3 YRa 
Pim = >—(Wo1To1 — Wim,out Tim) (1) 
Vim 
Rar, 
Pom = = ™ Wea,out a Wom, out) (12) 
om 


2.3. Thermal system 


The thermal system serves to reject the excess heat produced 
by losses in the stack. The major components of the circuit are a 
radiator for heat exchange with the ambient, and a fan to increase 
effectiveness of the heat by convection and radiation, and a water 
reservoir to store and to insulate the coolants thermally. Finally, 
a pump serves to supply the coolant to the heat source. The 
models for these components are based on the principle of heat 
transfer. Kroger [16] proposed a heat transfer coefficient, hrad, 
(kW m~? °C!) and pressure drop, pry (kPa) of a radiator as a 
function of the air mass flow rate Wai, (kg s7!): 


haa = —1.4495W2.. + 5.9045 Wair — 0.1157 (13) 


air 
Pr = (326.12 Wair — 75.396) + 101.325 (14) 
If the heat of the coolant is fully transferred to the radiator 


without any losses, the heat capacity of the coolant is identical 
with that of the radiator. Therefore, the outlet temperature of the 
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Fig. 5. Air and coolant controls. 


coolant can be expressed as a function of the radiator geometry, 
the temperature drop between the inlet and outlet, and the heat 
convection caused by the temperature difference between the 
ambient and the outgoing air temperature [16]: 


Frarea(Trad,c,in) Arad ) 
W.Cp, 


Trad,c,out = Trad,c,in —0.5 ( (15) 


where Frarea denotes the frontal area (m2) of the radiator, and 
Trad,c, in denotes the stack outlet temperature of the air. And the 
electric power for the fan can be calculated according to a ther- 
mal dynamic relationship between pressure drop and air flow 
rate [12]: 


Pran = (WaitC Pair Tams p 0/91) (16) 


Nelec "fan 


where Pfan denotes the electric power (W) and 7 does the effi- 
ciency of the fan and the fan motor. 


300 
250 
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100 


Stack current, A 
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The reservoir should be insulated thermally after the heat 
exchange occurs at the radiator by the convection. The variation 
of the heat in the reservoir is the sum of the heat that the coolant 
carries and the heat exchanging with the ambient. Accordingly, 
the outlet temperature at the end of the given time interval can 
be expressed by the following equation [12]: 


At 
MresCP yes 


F hA(Tyes,in = amb)) 


(W-Cp.(Ties,p = Tres, in) 


Ties, out = Tres,p ai 


(17) 


where Tes is the outlet temperature of a reservoir at the end 
of a time step (K), Tres,p the temperature of the reservoir at the 
previous of time step (K), At the time interval (s), Mres an equiv- 
alent mass of the coolant in the control volume (kg), hA the heat 
transfer of plumbing to ambient and Tres,in is the temperature of 
coolant into reservoir (K), respectively. 


O2 excess ratio 


15: T T 
19.0 195 20.0 205 21.0 215 22.0 
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Fig. 6. Comparison between sFF and sFB (a) load current; (b) oxygen excess ratio. 
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The coolant pump should deliver the coolants from the reser- 
voir and at the same time from the stack to radiator. If it is 
assumed that all heat energy generated in the stack is completely 
transferred to the coolant, then the mass flow rate of the coolant 
is expressed by its relationship with the heat source. 


Ò sou 


W. = 
° Cp. AT 


(18) 


2.4. Controls for the air and coolant flow rate 


Fig. 5 depicts a system configuration that includes two con- 
ventional PI controllers to control the flow rate of the air and the 
coolant. 

Pukrushpan et al. [17] proposed different control strategies 
for the air supply system to maintain the optimal oxygen excess 
ratio to prevent an oxygen starvation during abrupt changes of 
the load current. Hence, the oxygen excess ratio is defined as the 
oxygen supplied by the one consumed. Two control strategies, a 
static feed-back (sFB) and a static feed-forward (sFF) have been 
implemented and compared in Fig. 6. It turns out that the sFB 
demonstrates better performance with respect to the marginal 
limit of an oxygen starvation and a quick recovery behavior. 
Hence, the working temperature of the cell is assumed to be 
constant. 

When the load current varies abruptly, the temperature in the 
stack rises, particularly in the catalyst shown in Fig. 7. In addi- 
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Current, A 
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tion, Fig. 7 shows the temperature in the catalyst and coolant 
channel, where the temperature of the coolant is well controlled. 
According to the analyses conducted, the temperature in the 
catalysts on the cathode side is 3—8 °C higher than the average 
temperature of a stack dependent upon load current. The tem- 
perature in the catalysts is not measurable. Therefore, the actual 
temperature of the stack necessary for the controls is taken 
from the one at the outlet of the anode and the cathode side that 
has been averaged. The reference temperature of the stack is set 
to 77°C. 

It is observed that the oxygen excess ratio at the cell with 
varying temperature has increased, which results from different 
working temperature between the empirical mode and the model 
proposed. The temperature in the gas channel is low because of 
the cooling effect. Thus, the oxygen excess ratio shows a value 
of 2.07 with an offset. 

In addition, the oxygen excess ratio varies inversely as 
the current increases or decreases. At an increased current, 
more heat in the catalyst is generated and consequently the 
temperature in each of layers rises. In fact, the temperature 
of the coolant is the same as the reference value. Therefore, 
the temperature in the flow channel tends to rise as the 
load current increases and the partial pressure in the gas 
channel rises. Consequently, the mass of the air becomes 
less and the associated oxygen excess ratio gets lower. Like- 
wise, the oxygen excess ration gets higher as the current is 
decreases. 


yi T 
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Fig. 7. (a) Step current; (b) catalyst and coolant temperature; (c) oxygen access ratio. 
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3. DC/DC converter 


A typical configuration of a fuel cell powered system with a 
DC/DC converter is shown in Fig. 1. The PWM DC/DC con- 
verter can be described by using two state Eqs. (19) and (20), 
if the operation is limited to a conduction mode and no losses 
occur at the power conversion [19-21]. 


dupc 1. 1 ; VDC 
= AE 1—D 19 
dt Geer T. )iFC CRioaa (19) 
Ps Gees ) (20) 
= v v = —(v — v 
dt L FC DC L FC 1 


where vı denotes the average voltage at the lower switch during 
pulse width modulation (PWM) period, D the switching on-duty 
period, vrc the output voltage of the fuel cell system and vpe 
the output voltage of the converter, respectively. 

With no losses in the converter, the input power of the con- 
verter, Prc = ipcvpc, is identical to the output power of the 
converter, Poon = iconUpc. Then, the current at the output and 
the capacitor result is as follows: 


For the current loop, a proportional controller and integrator 
(PI) are used. The control variable for the loop is the current in the 
capacitor, which ultimately facilitates a complete control of the 
DC voltage. It is necessary to consider the dynamic of the load 
current, where the capacitor current is the difference between 
the output current and the load current drawn. Therefore, the 
load current is added, as shown in Fig. 5. The equations for the 
controls are given as follows: 


ss Kit \ g : 
v] == Rpt + y (cap ~ İcap) + VFC (24) 


, ki2 
ip = (Ho + =) (Ube — vpe) (25) 


The gains for the current loop have been particularly designed 
to increase the bandwidth of the closed-loop as high as possible 
to dynamically respond to any load changes. Conversely, the 
gains for the voltage loop are designed by considering the slow 
dynamics of the fuel cell system as a power source: 


UDC Gi+G2 
G(S)closed loop = —%_ = A (26) 
DC litoad=0 


where 


A = LCs* + KpE — D)s? + (Kp Kp + KRO — D)? + (Kp Kn + Kp2 Ki) — D)s + Ki Ki2( — D), 
Gi = KpK p — D)? + (1 — DXK pı Kin + Kp2Kia)s + Ka Kix — D), 


G2 = D(1 — D)K p8? + D0 — D)Kys. 


Pec UBCIFC 


icon = — = —— (21) 
UDC UDC 
leap = Icon — Hoad (22) 
; VECIFC 
Icon = ——— — Hoad (23) 
UDC 


The equations above present a second order linear system 
with two state variables, the DC output voltage vpc and the fuel 
cell current ipc. A block diagram for the controls of the PWM 
DC/DC converter is shown in Fig. 8. For the voltage control 
loop, a feed-forward controller is employed to increase con- 
trol dynamics on the voltage variation which directly feeds the 
reference voltage to the reference for the modulation. 


Feed-forward 
Controller 


DC/DC Converter 


4. Power control of the PEM fuel cell system with the 
DC/DC converter 


The objective of the power flow control is to accurately supply 
the power required from the load, thereby maximizing overall 
system efficiency and at the same time maintain the voltage 
level of the DC at the capacitor, which should not deteriorate 
the performance of the motor drive system. On the other hand, 
there are two energy sources available that can be combined to 
operate efficiently [22,23]. The total power that the system can 
supply is a sum of the energy stored in the DC capacitor and 
the power generated by the stack, which can be expressed as 
follows: 


Fig. 8. Control block diagram for the PWM DC/DC converter. 
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Table 2 

Simulation parameters for the DC/DC converter 

Switching frequency (kHz) 20 

Voc (V) 500 

C (uF) 1500 

Cin (WE) 350 

L (uH) 240 

Rioad (2) 6.25-62.5 


If a ratio of both power is defined as a new variable, s, then, 


Protal(t) = Prc(t) + Poc(t) = sPpc(t) + Poc(t) 
= (s+ 1)Poc() (27) 


The variable, s, determines a factor responsible for the power 
split between the fuel cell and the capacitor, which can be con- 
trolled according to the efficiency and the power being supplied 
by the capacitor. 

When the derivative of the DC voltage is negative, the capac- 
itor needs to be recharged from the stack. Otherwise, the DC 
capacitor is able to supply only a certain amount of power to the 
load. Then the reference power for the fuel cell stack is adjusted 
by a new value of the variable, s, which is determined by an 
efficiency map of the stack stored. Therefore a new reference 
for the stack current, Jf, is generated by a division of the DC 
power and the output voltage of the stack, which is expressed as 
follows: 


Pic = (sx(t) + s(x(t) — 1) Poc 


- UDC 

ce 0: ca <0 (28) 
(1+ s(t) Pp 
Bus DC 29) 
FC V, 
FC 

= 500 
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Fig. 9. DC voltage with the power control. 


IGBT forward voltage (V) 1 

IGBT off current tail time (s) 4x 10~° 
IGBT on (s) 0.017 
Diode on (Q) 0.015 
Diode forward voltage (V) 0.8 


Calculation step size (s) 


5. Integration and simulation 


Integration of this power system is conveyed on a real- 
time system in order to reduce the computational time. The 
fuel cell system is connected through a voltage source control 
block, where the current at the inductor is set to the reference 
current for the stack that controls the air flow rate. Subse- 
quently, the output voltage of the stack is adjusted according 
to the dynamics described. The output voltage and power of the 
stack at the operating power approximately 250 V. The mod- 
els developed are implemented with the software package of 
MATLAB/Simulink/Simpower. The parameters used for the 
simulation are summarized in Table 2. 
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Fig. 10. Dynamic response of the fuel cell power system at a step load with the 
power control. 
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Fig. 9 demonstrates the effects of the power control when 
the load current dynamically changes from 40 to 70 A. The yel- 
low areas present the contribution of the energy stored in the 
capacitors. When a load is applied, the voltage at the capaci- 
tor drops and the stack immediately reacts to refill the voltage 
gap by supplying the current. At the same time, the capacitor is 
being charged and the voltage at the capacitor increases. Then, 
the supply of the power from the stack is reduced to operate the 
stack more efficiently by the energy stored in the capacity. The 
switching mechanism repeats until the voltage reaches to the set 
steady state. 

Fig. 10 shows the simulated results for the current, voltage, 
oxygen excess ratio, temperature in the cathode catalysts as well 
as the gas flow channel of a fuel cell stack. When a step-like cur- 
rentis applied, the current and voltage of the stack follow the load 
changes with a time constant determined by the air supply sys- 
tem. When the new power controller is employed, the reference 
values for the stack current and voltage is continuously shifted 
to get an efficient operating point. Accordingly, the current and 
voltage varies dynamically and subsequently the oxygen excess 
ratio as well. However, it turns out that the temperature rise in 
the catalysts during the transition gets slightly lower than the one 
without power controller, which indicates reduced losses in the 
cell. At the same time, the temperature in the gas flow channel 
gets lower even at the increased load. 


6. Conclusion 


The paper presented has attempted to integrate component 
models for a PEM fuel cell power system that includes a stack, air 
supply system, thermal circuit and a PWM DC/DC converter as 
well as the associated control strategies. The integration reveals 
a variety of possibilities to study design parameters of a fuel cell 
powered system, particularly to consider a realistic characteristic 
of the fuel cell system along with an air supply system and the 
associated controls. A summary of new findings obtained is as 
follows: 


e the dynamic model for the stack currently used for controlling 
a power system can be improved by adding water balance in 
the membrane, diffusion effect in the GDL and temperature 
effects; 

e the oxygen excess ratio is inversely influenced by the load 
profile because of the change of the partial pressure in the gas 
flow channel; 

e the dynamics and efficiency can be improved when the energy 
stored in capacitor is utilized to continuously shift the oper- 
ating point of the stack to a high efficient area; 

e real time simulation with the fuel cell stack and the DC/DC 
converter enables a reduction of the computational time in two 
orders of magnitude that provides a new way to implement 
the automotive power system and analyze interacting effects 


between the electric load and the individual components of a 
cell. 


Future work will include: (1) integration of a bypass valve in 
the thermal circuit, (2) improvement of temperature control to 
effectively reduce temperature rise in the catalyst, (3) compen- 
sation of the temperature effect on the oxygen excess ratio and 
(4) advanced power management taking into account the fuel 
cell dynamics and efficiency. 
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